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P ARTICLE COUNT CORRECTION METHOD AND APPARATUS 

BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to counting of particles. The method is particularly 
useful in the counting of a sample of particles having a high density and variability in size. 
5 Background Art 

In many industrial and medical applications there is a need to determine the 
concentration of particles contained in a liquid sample solution. This sample solution 
generally contains one type of particle, which is of some minimum size and larger. These 
particles are counted by passing the sample solution fluid with the particles to be counted 

10 through an aperture, which has an electric current passing through it. When these particles 
go through the aperture an electrical pulse is generated whose peak voltage is proportional 
to the size of the particle. This analog particle voltage pulse is applied to a comparator 
circuit, which will generate a +5 volt level particle count pulse when the comparators 
threshold is exceeded. The particle count pulse signal is normally at 0 volts until the 

15 particle crosses the threshold and then goes to +5 volts for as long as the signal is above 
the threshold and returns to zero volts when the particle is no longer in the aperture 
sensing zone. This +5 volt particle count pulse then drives a counter, which is incremented 
for each particle count pulse on the leading or trailing edge. 

This counting technique works very well for lower concentration samples but is 

20 very inaccurate when high concentration samples are counted. There is a tendency to 
under count these higher concentration samples. This under counting is mainly caused by 
particle coincidence in the aperture. Coincidence occurs when more than one particle is 
sensed in the aperture at the same time, only one particle is counted. This results in 
undercounting for high concentration samples. 

25 Typically, to correct for this coincidence loss a coincidence correction algorithm is 

applied to the raw count to improve the count for high concentration samples. One of the 
most basic and effective statistically derived coincidence correction formulas is (Wynn 
and Hounslow - Coincidence correction for electrical-zone (COULTER COUNTER®) 
particle size analyzers, Powder Technology, 93, 163-175 (1997), see, in particular, p. 164, 

30 Formula 1. The entire contents of this reference are hereby incorporated by reference). 
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TC = True coincidence corrected Count in counts/s 



RC = Raw Count 
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ESV = Effective Sensing Volume in u.L 
FR = Flow Rate in uL/s 



Z = ESV/ FR in seconds 
The value of Z is adjusted until the result is linear. Parameter Z is a function of the 
aperture size, flow rate, and average particle size. 
m FC =TC * (DR I FR)*CF 



DR = Dilution Ratio 
CF = Calibration Factor 

This correction improves the count some but the loss at high concentrations is still 
too great and at higher concentrations the raw count will not increase with increasing 
15 concentration and may even recover lower counts. This coincidence correction formula is 
no longer effective when this happens. If the raw counts for these very high concentrations 
could be corrected then that would allow the coincidence correction algorithm to perform 
better. Another method is needed that will reduce the particle losses for higher 
concentration samples. 

20 The sample may be further diluted and rerun to reduce the loss and improve the 

raw count but this will result in reduced lab throughput as the sample has to be run again. 
The instruments dilution ratio may be increased to reduce the loss and improve the raw 
count but this will result in longer count times, which will impact the instruments 
throughput. It also degrades the count accuracy for low count normal samples. This would 

25 not be a preferred method to improve the count linearity. 

The optimum solution is to achieve an accurate and linear particle count for high 
and low concentration samples in the shortest amount of time. This is very difficult for 
some sample types such as white blood cells where the ratio of abnormal to normal cell 
concentrations is very high (450,000/5,000 or 90 to 1). The dilution is typically low (251 
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FC = Final Count in Counts/uL 
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to 1) to improve the throughput and accuracy of white normals. This low dilution makes 
the high concentration whites very difficult to count. 

Another method of count correction is well known in the industry (see U. S. Patent 
4,447,883 to Farrell et al, the entire contents of which are hereby incorporated by 
5 reference). This method calculates the average time between particles or total wait time of 
the counted particle stream and sets the corrected count equal to the inverse of the average 
wait time. The basic theoretical formula is that the frequency or counts per second is equal 
to the inverse of the average period of the particle pulse stream. 

This method assumes that the effects of the aperture size, particle size and flight 
10 time can be ignored. This method operates based on the wait time containing information 
about lost counts due to coincidence. If two particles go through the aperture too close 
together, then one count is lost. However, the two particles take up to twice as long to go 
through the aperture. This reduces the Total Wait Time and causes the true count to be 
increased. 

15 The Total Flight Time as measured in this method is corrupted with numerous 

coincidence events. Therefore the average flight time cannot be accurately determined 
from the Total Flight Time divided by the raw count. 

The Total Wait Time is an accurate measurement of the time that no particle is 
present. Therefore the average wait time derived by dividing the Total Wait Time by the 

20 raw count is a substantially reliable measurement. Also the Total Weight Time can be 
calculated by measuring the total flight time and subtracting this from the total count time. 

One key parameter is not handled very well in current methods including both 
methods mentioned above. This is the consideration of the effects that particle size have 
on proper coincidence correction. Both methods above would not work properly if 

25 different size particles were run through the aperture. The coincidence correction method 
uses a fixed average particle size to determine the effective sensing volume. If the particle 
size changes from sample to sample then the formula will no longer work properly. The 
wait time count correction method does not include the true Total Flight Time in the 
formula. Even though it measures the Total Flight Time it only uses this information to 

30 calculate the Total Wait Time. The Total Flight Time as measured is corrupted by 
coincidence and is not a good measurement of the true flight time. This method also will 
not perform well with different size particles. This is important for sample types such as 
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white blood cells. White blood cells typically have a wide range of particle sizes. Clearly, 
the effects of particle size must be included in the correction methods to make them work 
properly at high concentrations with multiple size particles. Most count systems also 
generate a size histogram for the particles being counted. The average size of the 
5 population counted can be determined from the histogram. 

Various methods for coincidence correction are described in Wales and Wilson, 
The Review of Scientific Instruments, Vol. 32, No. 10, p. 1132-1136 (1961), Princen and 
Kwolek, The Review of Scientific Instruments, Vol. 36, No. 5, p. 646-653 (1965), U. S. 
Patents 4,580,093 to Feier et al, 5,247,461 to Berg et al, and 5,452,237 to Jones, Jr. 
10 Most count correction methods rely on only one technique to perform the 

correction. It would be better to utilize more than one technique to allow for more gentle 
correction from each one with a more linear (accurate) final result. 

Summary of the Invention 

This present invention relates to an improved particle count correction method with 
15 better particle count linearity (better accuracy) for samples with wide variation in particle 
sizes and wide variations in concentrations (high concentrations often produce lower 
accuracy). 

This is accomplished by providing a method for counting particles, comprising the 
steps of successively passing multiple particles through a particle sensing zone; 

20 introducing a first signal into said particle sensing zone for a period of time; measuring a 
second signal emanating from said particle sensing zone, said second signal being caused 
by modulation of said first signal by said particles passing through said particle sensing 
zone; generating raw data using said second signal, said raw data correlating to a raw 
count of particles passing through said chamber, a wait time count and a size of each 

25 particle; and processing said raw data by using a true average flight time and a true 
average wait time to obtain a corrected count of particles. 

In another embodiment, the invention is directed to an apparatus for counting 
particles in a sample, comprising one or more particle sensors, each sensor having a 
sensing zone; a particle delivery unit for delivering particles to at least one of said particle 

30 sensing zones, said particles passing through at least one sensing zone; a particle 
measuring unit for determining the size of particles passing through at least one of said 
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particle sensing zones, said sensor generating a particle size signal, and for determining 
the number of particles that pass through at least one of said particle sensing zones in a 
given period of time, said particle sensor generating a particle number signal; a device for 
calculating the average flight time of said particles in said sample based on said particle 
5 size signal and said particle number signal; and a correcting unit for correcting an apparent 
particle count to an adjusted particle count by adding a true average flight time to a true 
average wait time to obtain a corrected count of particles. 

In still another embodiment, the invention is directed to an apparatus for counting 
particles, comprising a chamber having an inlet, an outlet and a particle sensing zone 

10 located between said inlet and said outlet; a pump for passing a fluid containing particles 
into said inlet, through said particle sensing zone and out of said outlet; an electric source 
arranged to pass an electric current through said particle sensing zone; an electric current 
detector for measuring electric current as particles pass through said particle sensing zone, 
said detector generating raw data indicative of the number of particles passing through 

15 said particle sensing zone and indicative of the size of particles passing through said 
particle sensing zone; and a program for processing raw data from said detector, said 
program having the capability add true average flight time to average wait time to give a 
true average period value. 

In yet another embodiment, the present invention improves the wait time count 

20 correction method by adding the true average flight time to the average wait time to give a 
true average period value which will yield a more accurate count for samples with a wide 
range of cell sizes such as white blood cells. This is more accurately an average period 
count correction method. The true average flight time is derived from the mean value of 
the size histogram of the sample population and a formula which relates the particle flight 

25 time to mean size. This formula is empirically derived from the counting instrument and is 
a function of aperture size, particle size, and flow rate. 

This method also uses both the Average Period Count correction method and an 
enhanced coincidence correction formula to correct the raw count. Adjustment factors 
FTCF, BaseZ, and CCRFT are used in their respective methods to adjust the amount of 

30 correction each one has on the count. 

Average Period correction method formula: 
APC = Average Period Count 
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CT = Count Time in seconds 

RC = Raw Count for one CT or the average from multiple CT's. 
WTC = Wait Time Count for one CT or the average from multiple CT's. 
CTClk = Count Time Clock in counts per second 
5 AFT = Average Flight Time in microseconds, 

(Count conversion 10 is in units of counts per microsecond) 
FTCF = Flight Time Correction Factor 



WTC + FTCF * (RC * 10 * AFT ) 
(3 ) CTClk *CT 

Enhanced Coincidence Correction formula: 
10 TC = True Count 

APC = Average Period Count 

BaseZ = Main factor for control of coincidence correction 

CCRFT = Coincidence Correction Flight Time factor. 

Z = BaseZ + (AFT * CCRFT) (includes correction for particle size) 



TC = — ln( l - ) 

(4) Z y \-{Z*APCY 



Brief Description of the Drawings 

Fig. 1 depicts the various operational units of a preferred embodiment of the 

device; 

Fig. 2A is a pulse diagram showing the interaction of different size particles and 
coincidence in the aperture waveforms and identifies the corresponding flight time pulse 
width and wait time, as well as the corresponding count and wait time count; 

Fig. 2B depicts the Digital Count pulse. 

Fig. 2C depicts the Clock pulse. 

Fig. 2D depicts the Wait Time Count Clock pulse. 
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Fig. 3 is a graph that shows the improved linearity of one embodiment of the 
inventive method when compared to another conventional method using a readily 
available WBC linearity control material; 

Fig. 4 is a graph that illustrates the average Raw Count (ARC) values and the 
5 results obtained after applying each method using the same data as Figure 3. APC shows 
the effect that the average period method has on ARC. TC shows the results of applying 
the enhanced coincidence correction method to APC; 

Fig. 5 shows the correlation between the method of the invention and a reference 
method for abnormally high WBC blood samples. The reference method values were 
10 obtained by pre-diluting the samples and running them on a Beckman Coulter, Inc. 
GEN*S™. The values using the present invention were obtained from a Beckman 
Coulter, Inc. LH-750™ instrument without dilution. 

Figs. 6-9 are WBC histograms for four different samples from data in Fig. 5 having 
an abnormally high WBC count. 

15 D etailed D es cription of the Invention 

Definitions 

In the present application, unless the context indicates otherwise, the following 
terms and abbreviations have the following meanings. 

Average Period Method: A method for correcting particle counts where high 
20 levels of coincidence loss are present. It generates the corrected particle count from the 
inverse of the derived average period of the particles counted. It will automatically correct 
the count when different size particles are run in the instrument. This method performs 
well until the particle concentration becomes very high. 

Enhanced Coincidence Correction Method: A method for correcting particle 
25 counts where high levels of coincidence loss are present. It generates the corrected particle 
count from a standard coincidence correction algorithm that has been enhanced to include 
the average flight time to allow adjustment of the amount of correction required based on 
particle size. It will automatically correct the count when different size particles are run in 
the instrument. This method performs well by itself on the raw count until the 
30 concentration becomes very high. 



-8- 

Average Period Count (APC): The corrected particle count derived from the 
average period correction method. This uses the raw count, the wait time count, and the 
size histogram to correct the raw count. 

True Coincidence Corrected Count (TC): The corrected particle count derived 
5 from the enhanced coincidence correction method. To improve linearity and accuracy, it is 
applied to the Average Period Count. 

Final Count: The true count is multiplied by the scaling factor (DR/FR) and 
calibration factors to yield the final particle count result in counts per uL. 

Correction Factors ( FTCF, BaseZ, CCRFT ) : These correction factors are 
10 empirically determined in parallel by using various size latex bead solutions at various 
concentration levels as well as blood control material. Only one size of bead is used per 
test solution. They are adjusted until the count linearity is optimized for the various size 
test particles. For any instrument there are numerous sets of proportional values, which 
yield acceptable results. Different instruments will require different factors depending on 
15 the aperture size, particle size, and flow rate. Some example of correction factors include: 
Flight Time Correction Factor (FTCF) : This factor adjusts the amount 
of correction performed by the average period method. 

BaseZ : This factor adjusts the amount of correction performed by the 
coincidence correction algorithm. 
20 Coincidence Correction factor from the Flight Time (CCRFT) : This 

factor adjusts the amount of correction of the coincidence correction algorithm 
using the derived average flight time. 

Count Time (CT): The time that particles are counted in the aperture. This method 
has a preferred count time of one second that is repeated from 4 to 10 times depending on 
25 whether the histogram sizing is complete. However, other suitable count times well 
known in the art may also be used. Having information about counts per second is very 
useful in debugging problems. Also, it is possible to eliminate counts that are significantly 
different than the rest. This eliminates intermittent problems due to bubbles or clogs. 
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Raw Count (RC or ARC): The total count in a counter that is incremented by the 
rising or falling edge of the count pulse for exactly 1 count time. The raw count is 
depicted in Fig. 2B. This counter is buffered and reset at the end of each second. The 
system then reads each of the buffered values. If the count time is repeated multiple times 
5 then this value can be the average of the multiple counts. 

Count Time Clock (CTClk): The clock used to increment the wait time counter 
when no particle is present in the aperture sensing zone. This is in units of counts per 
second. For example, a 10MHz clock CTClk = 10,000,000 counts/second. An example of 
a graphical depiction of the 10MHz clock pulse is shown in Fig. 2C. It is used to convert 
10 the wait time count into seconds. 

Wait Time: The time spent waiting for the next particle after the current particle 
has exited the aperture sensing zone. This is the time between two particles when no 
particle is sensed in the aperture and the count signal is low, as is shown in Fig. 2A. 

Wait Time Count (WTC): The total count from a counter that is enabled to 
15 increment the count from a fixed frequency clock when no particle is present in the 
aperture during the count time. If the count time is repeated multiple times then this value 
can be the average of the multiple counts. A graphical depiction of the Wait Time Count 
with respect to the digital count pulse and Count Time Clock is shown in Figs. 2B-2D. 

Total Wait Time (TWT): The total time that no particle of sufficient size is in the 
20 apertures effective sensing volume during the count time. This is the sum of all the spaces 
between the analog pulses that cross the threshold. TWT = WTC / CTClk. This is a 
measurement of the total time between the particles as it is not corrupted at high 
concentrations like the flight time. 

Flight Time (Pulse Width) (FT): The time that a particle is sensed in the 
25 apertures effective sensing volume such that the resulting analog pulse is above the count 
threshold. As is shown in Fig 2A, this is the time that a particle of sufficient size is sensed 
as it flies through the aperture. This is equivalent to the pulse width time of the analog 
pulse as measured at the points where the analog pulse goes above and then below the 
count threshold. 

30 Flight Time Count (FTC): The total count from a counter that is enabled to 

increment the count from a fixed frequency clock when a particle is present in the aperture 
during the count time. 
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Total Flight Time (TFT): The total time that a particle of sufficient size is in the 
apertures effective sensing volume during the count time. This is the sum of all the analog 
pulse widths that cross the threshold. The typical formula is TFT1 = FTC/(CTClk*CT). 
This measurement is corrupted by coincidence events at high concentrations and is not an 
5 adequate measurement of the true total flight time of the particles. This new invention 
method uses the Average Flight Time and the Raw Count to provide a more accurate 
estimate of the true total flight time. TFT2 = (RC * 10 * AFT)/(CTClk*CT). 

Average Channel (AC): This is the channel number (0 to 255) derived from the 
256 channel size histogram. This channel number is representative of the average size of 
10 the particles that were counted. 

Average Flight Time (AFT): This is the estimated true average flight time or 
pulse width of the particles that were counted. This is the average time that a particle pulse 
was above the count threshold. This is derived from the average channel and the 
empirically derived formula that converts the channel number to pulse width or flight 
1 5 time. It does not matter if the particles in a sample are of one size or have a wide size 
variation the average flight time value will linearize the count. 

Z: This is the main adjustment factor used in the coincidence correction formula 
used in the method according to this invention. This is essentially equivalent to the 
average time that a particle is present in the aperture sensing zone. It is typically in units of 
20 micro seconds. This is determined by the ESV or effective sensing volume and the FR or 
Flow rate. The formula is Z(s) = ESV(uL) / FR(uL/s). 

Effective Sensing Volume (ESV): This is a volume in micro liters that consists of 
the effective aperture volume plus the increase in the effective volume due to the particle 
size. Because the particle is sensed before it is actually in the aperture the effective 
25 aperture volume looks larger. Larger particles will be in the sensing zone longer than 
smaller particles. This part of the volume is variable depending on the size of the particles. 
Typically, the average size of the particle is used. The enhanced method described does 
include correction for particle size variation. 

Particle samples 

30 The present invention has particular application to counting particles in samples 

having high density (or a potential for significant density variability) and/or high particle 
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size variability. For example, the invention is particularly useful for samples that have a 
particle concentration so high that the average time between particles is less than the flight 
time. Cell density as it flows through the aperture is ((450,000counts / uL) / 251)/ 
1.12(Cal factor) = 1600 cells per microliter. White cells vary in size from 30 to 450 
5 femtoliters. In the LH-750™ this corresponds to flight times of from 15 to 30 
microseconds, respectively. The WBC flow rate is 40 microliters per second. The average 
particle count rate for a 450,000 WBC count blood is Count Rate = ((WBC * FR) / DR) / 
Cal Factor = ((450000counts/uL * 40uL/s) / 251) / 1.12 = 64,029. This corresponds to an 
average cell period if the cells went through the aperture evenly spaced of 1/ 64,029 = 15.6 
~;i 10 microseconds. At this high concentration it is clear that if the cells went through evenly 
;! spaced that only one cell would be sensed as they are too close together. The pulse width 

; !: | of the white cell is equal or greater than the average spacing between the cells. It is only 

1: because the cells are randomly spaced that a count can still be recovered. When both high 

W density and high particle size variability are present in a sample, current methods can often 

15 render inaccurate particle number counts. However, the invention also works with 
J;- samples having a substantially homogenous particle size distribution. 

M In a preferred embodiment, the particles to be counted are blood cells such as white 

rT blood cells and red blood cells. However, the present invention is applicable to counting 

many different types of particles, especially where particle size variability and particle 
20 density variability are expected. Examples of other types of particles that may be counted 
include latex particles and other biological or industrial particles. The preferred 
embodiment of the invention applies to counting single types of particles or cells of a 
certain minimum size and larger. 

Counting platelets is especially difficult in that they usually have larger red cells 
present at the same time, which can significantly reduce the PLT count. The red cells 
cause a coincidence loss in the PLT counts when a PLT goes through the aperture with a 
red cell. In order to apply this invention to platelet counts when red cells are counted in the 
same bath, formula (4) needs to be modified to reflect the coincidence loss due to the 
presence of the red cells. 

White blood cells are usually more variable in size than red blood cells and 
25 platelets, and therefore, the present invention has particular applicability in counting white 
blood cells. For example, in a "normal" blood sample (a sample from an individual that 
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does not have an infection), white blood cells often vary in size by more than 50%, and 
can vary by more than 100% (with the smaller blood cells representing the baseline size 
for measuring variability). In abnormal bloods the white cells size can vary greatly 
depending on the nature of the illness. Some samples have small white cells as in Fig 9 
5 and some have large cells as in Fig 8. See also Figures 6 and 7. Thus, the present invention 
allows more accurate counting of samples that have a wide variation (50fL to 400 fL) in 
the size of the white cells. 

In addition, the number of white blood cells (density) can vary greatly from a 
"normal" (uninfected) patient to an "abnormal" (infected) patient. For example, the 

10 number of white blood cells in an infected patient can be 90 times greater than the number 
of white blood cells in a normal patient. Thus, the present invention allows more accurate 
counting of samples that are expected to have concentration differences of 2 fold or 
greater, preferably 10 fold or greater, more preferably 75 fold or greater, and even 100 
fold or greater. Thus, the particle counter can accurately count white blood cells in a low 

1 5 density sample (from a normal individual) and in a high density sample (from an infected 
individual) under the same testing conditions. This prevents having to rerun a sample 
through a machine (after dilution) if abnormally high numbers of cells are detected in 
order to obtain a more accurate reading. 
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Preferred Embodiment 

As shown in Figure 1, a sensor in a particle counting device 2 makes 
measurements relating to the size of particles in a particle sensing zone and the number of 
particles that pass through the particle sensing zone in a given period of time. The sensor 
may be in a particle counting device 2 such as a COULTER COUNTER® or other 
5 suitable particle counting device well known in the art. In such a device, particles 
suspended in an electrically conductive medium (usually an aqueous salt solution) are 
measured by passing the particle suspension through a small orifice or aperture through 
which an electric current is flowing. This is the particle sensing zone. A particle passing 
through the particle sensing zone produces a decrease in conductivity; the sensor converts 

10 the conductivity changes to a voltage pulse (particle pulse signal), which can be counted 
(see Figs. 2A and 2B). The magnitude (amplitude) of the pulse is nearly proportional to 
the volume of the particle. If two particles pass through the particle zone at the same time 
(or too close together), the two particles can be counted as a single particle (See Fig. 2A). 
This is known as "coincidence loss". 

15 Each COULTER COUNTER or other particle counting device may have one or 

more apertures and a sensor corresponding to each aperture. 

The sensor in the particle counting device 2 is connected to a preamp 4, which 
receives the output particle pulse signal generated by the sensor. The preamp 4 amplifies 
the signal and generates an analog voltage signal Vsig 6. The preamp 4 is also connected 

20 to a comparator 8, peak detector 26 and control and editor circuit 30. 

The comparator 8 is connected to the preamp 4, AND gate 20, and particle counter 
12. The comparator 8 receives Vsig 6 from preamp 4 and compares Vsig 6 to a threshold 
Vth (shown in Figure 2A) and generates digital signal 10 representing the portions of Vsig 
above the threshold Vth. The digital signal 10 of the comparator 8 is received by the 

25 particle counter 12 and the AND gate 20. 

The particle counter 12 is connected to comparator 18 and average raw count 
generator 14. The particle counter 12 generates a raw count that is incremented by the 
rising or falling edge of each pulse in signal 10. 

The average raw count generator 14 receives the output raw count from the particle 

30 counter 14 and generates an average raw count ARC. 
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A clock 18, which is preferably a 10 MHz clock, is connected to the AND gate 20. 
The AND gate receives signal 10 from comparator 8 and signals from the clock 18. The 
AND gate 20 outputs a signal from the clock 18 when digital signal 10 allows the clock 
signal to pass through the AND gate 20. 
5 The wait time counter 22 is connected to AND gate 20 and the average wait time 

count generator 24, and when the clock signal is output from the AND gate, the wait time 
counter 22 increments to generate a wait time count. 

The average wait time counter 24 is connected to the wait time counter 22 and the 
average period count generator 16. The average wait time count generator 24 receives the 
10 wait time count from the wait time counter 22 and generates an average wait time count 
AWTC. 

The peak detector 26 is connected to preamp 4, control and editor circuit 30 and 
analog to digital converter 28. The peck detector receives Vsig 6 from preamp 4 and 
outputs the amplitude of a peak value of the pulses Vpk to the analog to digital converter 
15 28. Vsig and Vpk are also outputted to the control and editor circuit 30. 

The control and editor circuit 30 is connected to the peak detector 26 and buffer 
32. The control editor circuit 30 receives and compares Vpk and Vsig to edit the pulses by 
suitable techniques well known to those of ordinary skill in the art. The control and editor 
signal generates a write pulse for good cells to store only the size data of good cells in data 
20 buffer 32. 

Digital converter 28 is connected to buffer 32 and receives the Vpk data and 
generates digital size data, which is stored in data buffer 32. 

Data buffer 32 is also connected to the channelizer 34, which may read the data 
stored in buffer 32, such as the digital size data of good cells, to generate a channelized 
25 size histogram. Examples of the channelized size histograms are provided in Figures 6-9. 

The channelizer is also connected to smoother 36, which receives the channelized 
size histogram from the channelizer 34 and smooths the histogram. 

A histogram mean and mode generator 38 is connected to the smoother 36 and 
average channel generator 40. The histogram mean and mode generator 38 receives the 
30 smoothed histogram from the smoother 38 and generates an estimate of the average 
population size. 
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The average channel generator receives the estimate of the average population size 
from the histogram mean and mode generator 38 and outputs an average channel size AC. 

The average flight time generator 42 is connected to the average channel generator 
40, the Z factor generator 46 and the corrected average flight time generator 44. The 
5 average channel size AC is received as an input to the average flight time generator 42 and 
the average flight time generator 42 converts the average channel size into an average 
flight time per pulse AFT. 

The average flight time AFT is then received by both the corrected average flight 
time generator 44 and the Z factor generator 46. The corrected average flight time 
10 generator 44 multiplies the average flight time AFT by the flight time correction factor 
stored therein and generates the corrected average flight time CAFT. 

The average period count generator 16 is connected to the corrected average flight 
time generator 44, average raw count generator 14 and the average wait count generator 24 
and uses the outputs of all three generators 44, 14 and 24 to calculate an average period 
15 count APC. 

The Z factor generator is connected to the average flight time generator 42 and the 
coincidence corrected count generator 48 and receives the average flight time AFT from 
the average flight time generator 42 and also receives inputs of the coincidence correction 
factor from the flight time CCRFT and BaseZ factor. Using these inputs, the Z factor 
20 generator calculates the Z factor and outputs the Z factor to the coincidence corrected 
count generator 48. 

The coincidence corrected count generator 48 is connected to the average period 
count generator 16 as well as the Z factor generator 46, and receives the average period 
count APC output by the average period count generator 16 and the Z factor output by the 
25 Z factor generator 46. Using the average period count APC and the Z factor, the 
coincidence corrected count generator 48 performs coincident correction calculations on 
the data and outputs a true count TC. 

The true count is received by multiplier 50 which is connected to the coincidence 
corrected count generator 48. Multiplier 50 converts and scales the data to generate a final 
30 count 52 of the counted particles such as white blood cells in the sample. 

As discussed above, the particle sensor maybe part of a COULTER COUNTER or 
other suitable particle counting device 2 well known in the art. The preamp 4, comparator 
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8, particle counter 12, clock 18, AND gate 20, peak detector 26, analog to digital converter 
28, and control and editor circuit 30 also may be part of the COULTER COUNTER or 
other suitable particle counting device as hardware components. The wait time counter 22 
may be a hardware or software component operated by microprocessor within the 
COULTER COUNTER or other suitable particle counting device. The average raw count 
generator 14, average wait time count generator 24, average period count generator 16, 
channelizer 34, smoother 36, histogram mean generator 38, average channel generator 40, 
average flight time generator 42, corrected average flight time generator 44, Z factor 
generator 46, coincidence corrected count generator 48, and multiplier 50 are software 
applications which are embodied and performed in a computer system, computer or 
microprocessor associated with the COULTER COUNTER or other suitable particle 
counting device. 

It is noted that the invention may alternatively have various embodiments as would 
be known to one of ordinary skill in the art. For example, instead of hardware, software 
applications could be utilized for each of the comparator 8, particle counter 12, peak 
detector 26, control and editor circuit 30, and wait time counter 22. Further, the software 
applications may be embodied in one or more computer processors or microprocessors. 

The following represents an example of how a preferred method of the operation of 
the invention can be applied to the WBC parameter in an LH-750 with three apertures. It 
can be applied just as well to other parameters, i.e., RBC, any number of apertures, and 
other suitable particle containing devices. 

The method is described with respect to one aperture of three apertures in the LH- 
750. However, similar processing is performed for each aperture. The data from each 
aperture is later combined by multiplier 50 as discussed below to generate the final count. 

Step 1: Get Raw Count and Wait Time Count Data 

As shown in Fig. 1, a sensor in a particle counting device 2 generates waveform 
signals (particle pulse signals) corresponding to particles passing through an aperture in 
the sensor (see also Figs. 2A-2D). The voltage waveform signal represents the particle 
passing through the aperture as a function of time. The voltage is affected by the size of 
the particles. The waveform signals from the sensor are fed through a preamp 4 to amplify 
the signal up to about 10 volts. An amplified analog voltage signal 6 (Vsig) is sent from 
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the preamp 4 to the comparator 8 where the analog signal is converted to a digital signal 
10. All signals that fall below a certain voltage threshold (Vth), as shown in Fig. 2A, are 
filtered out in the comparator 8. Signals 10 from the comparator 8 are sent to a particle 
counter 12 to generate a raw count where particles are counted for one second of count 
5 time. The average raw count generator 14 uses the raw count from the particle counter 12 
to generate an average RC. The average RC is calculated by averaging multiple raw 
counts (such as 4-10 raw counts). For example, four 1 -second raw counts or ten 1 -second 
raw counts can be averaged to generate an average RC. 

The average RC is then used as one of 3 inputs to generate an average period count 

10 (APC). Note: Any individual one second values may be excluded based on being non- 
uniform as compared to the other values. 

Digital signals 10 from comparator 8 are also used to enable the count time clock 
18 (i.e., 10MHz constant frequency clock) to pass through AND gate 20 and increment the 
wait time counter 22 when digital signal 10 is low or zero volts indicating that no particle 

15 of sufficient size is in the sensing zone. The value in the counter at the end of the count 
time is the Wait Time Count. The average wait time count generator 24 uses the wait time 
count from the wait time counter 22 to generate an average WTC. The average WTC is 
calculated by averaging multiple wait time counts (such as 4-10 wait time counts). For 
example, four 1 -second wait time counts or ten 1 -second wait time counts can be averaged 

20 to generate an average WTC. The average WTC is then used by the average period count 
generator 16 as one of 3 inputs to generate an average period count (APC). Note: Any 
individual one second values may be excluded based on being non-uniform as compared 
to the other values. 

This is different from the previously referenced U. S. Patent 4,447,883 at least in 
25 that the time between particles is directly measured instead of deriving the time by 
subtracting the measured total flight time from the total count time. 

Step 2: Correct the Raw Count value using the Average Period method. 

The average flight time (AFT) or average pulse width cannot be derived from the 
measured Total Flight Time and the raw count like the wait time can. The measured total 
30 flight time is corrupted because of the coincident particles that go through the aperture at 
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the same time. In order to calculate the average flight time the average cell size must first 
be estimated. 

In order to accomplish this, the analog signal 6 is sent to a peak detector 26 which 
stores the peak value (Vpk) or amplitude of any pulse. The amplitude of an analog pulse 
5 is proportional to particle size. The Vpk is sent to an analog to digital converter 28 where 
a digital signal is generated which is proportional to the analog Vpk value. The analog 
signal 6 and Vpk are also sent to a control and editor circuit 30 where real time signals are 
compared with peak signals to perform editing of the pulses. This editing process is well 
known in the industry and will eliminate those particles that have abnormal analog 

10 waveforms. See, for example, U.S. Patent 3,710,264 incorporated herein by reference. 
This helps make the size histogram a true representation of the sample population that was 
counted. A write data signal is only pulsed for good cells and therefore bad cells are 
rejected and not written to the data buffer 32. The digital size value from converter 28 is 
stored in the data buffer 32 when the write data signal from circuit 30 is pulsed. Data 

15 buffer 32 stores the size information of each good cell as determined by the control and 
editor circuit 30. 

The channelizer 34, smoother 36, histogram mean and mode generator 40, and 
average channel generator 40 operate using methods well known to those of ordinary skill 
in the art. 

20 The stored digital size data from the buffer 32 is read and channelized into a size 

histogram by channelizer 34. The channelizer 34 scales the digital size value and uses the 
scaled digital size to select and increment one of the bytes or channels in a 256 byte array. 
This generates the size or frequency histogram. Examples of size histograms are depicted 
in Figs. 6-9. 

25 The average mean size of all particles that went through the sensing zone can be 

calculated from the channelized histogram of all the particles as indicated at 34. Smoother 
36 smooths the size histogram using any suitable standard smoothing algorithm. 

Mean and mode values are generated from the smoothed histogram by the 
histogram mean and mode generator 38 using equations 5 and 6 below. However, any 

30 method known to those of ordinary skill in the art which yields a suitable estimate of the 
true population size can be used. In this example, the average of the histogram mode and 
mean channel is used. The average channel generator 40 uses equation 7 to generate the 
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average channel. This calculation corrects for the distortion in the histogram due to 
coincidence and signal conditioning. The goal is to get the best estimate possible of the 
true average size of the particles that are counted. 

(5) AMDC = Average Mode Channel = (Apl MaxCh(0..255) + 
Ap2MaxCh(0..255) + Ap3MaxCh(0..255) )/3 

(6) AMNC = Average Mean Channel = (AplMeanCh(0..255) + 
Ap2MeanCh(0..255) + Ap3MeanCh(0..255) )/3 

(7) AC = Average Channel = (AMDC + AMNC)/2 

The average flight time generator 42 converts the Average Channel to the Average 
Flight Time per pulse using an empirically derived formula. The relationship between 
particle size and particle pulse width should be characterized for a particular aperture size 
and flow rate. Various size particles are run through the instrument and the relationship 
between the particle size and pulse width is characterized. These values would be 
instrument related as each instrument may have different bandwidths and throughput 
which would affect the pulse width time. Each system may have a different characteristic 
relationship between the particle size and the Pulse width time. 

Various formulas can be used for Average Channel AC to Average Flight Time 
AFT conversion. Two exemplary formulas are: 

(8a) AFT = 15.47 + (AC/ 13.7) and 

(8b) AFT = 7.3983*Ln(AC) - 10.431 

The corrected average flight time generator 44 multiplies the Average Flight Time 
by the FTCF correction factor to yield the CAFT or Average Corrected Flight Time. This 
calculation is exemplified by equation 9. 

FTCF = Flight Time Correction Factor - 0.59 ( This was empirically determined 
with latex beads using the Beckman Coulter LH-750. This allows the amount of 
the correction to be adjusted. However, other suitable Flight Time Correction 
Factors may be utilized.) 

(9) CAFT = Corrected Average Flight Time = FTCF * AFT = 0.59 * AFT 
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The CAFT generated by the Average Flight Time Generator is one of three inputs that are 
used to calculate the APC value by Average Period Count Generator 16. 

The Average Period Count Generator 16 generates the APC values using the 
following formulas. The Average Period Count general formula is 
5 APC = ARC (Average Raw Count from 14) / APT( Average Period Time). 

The Average Period Time formula is 

APT = TWT (Total Wait Time) + TFT2 

where 

TWT = Total Wait Time = WTC( from 24) / (CTClk * CT) 

10 and 

TFT2 = Total Flight Time = (CAFT (from 44)* 10 * RC (from 14) / (CTClk*CT). 
Since the CAFT is in microseconds, it must be multiplied by 1 0 to get units of 1 00 ns 
which the wait time count is in. Alternately the APT formula can be written as 

APT = (WTC + (CAFT * 10 * RQ) / (CTClk*CT). 
1 5 The Final Formula that is used by the Average Period Count Generator 16 is: 

(10) APC = RC / ( (WTC + (CAFT * 10 * RQ) / (CTClk*CT) ) 
Step 3: Coincidence Correct the Average Period corrected Count. 

Further correction of the data is obtained by using APC and Z and performed by 
the Coincidence Corrected Count Generator 48. 
20 The Z factor is the main adjustment factor and is generated at by the Z factor 

generator 46 using equation 11. 

(11) Z = BaseZ + (AFT * CCRFT) = 34.35* 10 A -6 + ( AFT * -0.3*10 A -6), where 

BaseZ = Z base value = 34.35 *10 A -6 (empirically determined with latex beads) 
AFT = Average Flight Time from step 2 (42 in Fig. 1) 
25 CCRFT = Coincidence correction factor from flight time = -0.3 *10 A -6 

(empirically determined with latex beads) 



Although the Z base value and Coincidence Correction factors are empirically 
determined, these factors are examples of suitable correction factors, and other suitable 
factors may be used. 
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The coincidence corrected count generator 48 generates a True Count TC using 
equation 12. 



(12) Z 1~Z*APC , wherein 

TC = True coincidence corrected Count, and 
5 APC = Average Period Count 

Multiplier 50 receives the True Count TC from the Coincidence Corrected Count 
Generator 48. Although Fig. 1 indicates only one line 54 from the Coincidence Corrected 
Count Generator 48 to Multiplier 50, the Multiplier receives a True Count TC 
corresponding to each aperture. In the embodiment described, there are three apertures. 
10 The Multiplier 50 receives three True Counts. In other embodiments, the number of True 
Counts TC received by the multiplier could be more or less than three. 



Step 4: Scaling, Matching, Averaging and Final Count 

The Multiplier 50 performs calculations on the True Counts TC from each aperture 
15 to generate the final count 52. The calculations include scaling, matching and averaging 
the data by methods well known to those of ordinary skill in the art. 

The data is converted back to count/ul by multiplier 50, so that it is in a form that 
is convenient or usable. Other conversions to other forms may also be performed as 
needed. Values from multiple channels (if present) are averaged. A final count 
20 representative of the total number of white blood cells in the sample is obtained as shown 
at final count 52. 

The Average Period count correction method is more accurate than previous 
methods. It takes into consideration the effects of particle size in the correction formula so 
that it is linear, as is shown in Fig. 3 (which corresponds to accuracy), no matter what size 

25 of particle is run in the system. By using this method with an enhanced coincidence 
correction formula, the count linearity can be maintained even at very high sample 
concentrations such as are found in high 450K white blood samples. By allowing the 
concentration to be this high, the throughput of the instrument is kept high. Throughput 
and accuracy can be achieved with current technology which is adapted to use the above- 

30 described correction method. 
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The invention is not limited to the above embodiments. For example, in the above 
embodiment, averaging of multiple count times occurs on the raw data in the Average 
Raw Data Count Generator 14 and Average Wait Time Count Generator 24. These 
averaging steps could also be performed at the Average Period Count Generator 16 or 
5 after the Coincidence Corrected Count Generator 48 and even at the final count 52. Four 
to 10 Average Period Count APC Values could be calculated from the raw data and then 
averaged; four to 10 true count TC values could be calculated and averaged; or, four to 10 
final count values could be calculated and then averaged. 

The invention being thus described, it will be obvious that the same may be varied 
10 in many ways. Such variations are not to be regarded as a departure from the spirit and 
scope of the invention, and all such modifications as would be obvious to one skilled in 
the art are intended to be included within the scope of the following claims. 



